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We investigate the effect of unconventionai density wave (UDW) condensate on an Anderson 
impurity using large-iV technique at T = 0. fn accordance with previous treatments of a Kondo 
impurity in pseudogap phases, we find that Kondo effect occurs oniy in a certain range of param- 
eters. The f-eiectron density of states reflects the influence of UDW at iow energies and around 
the maximum of the density wave gap. The static spin susceptibiiity diverges at the criticai cou- 
pling, indicating the transition from strong to weak coupling. In the dynamic spin susceptibility 
an additional peak appears showing the presence the UDW gap. Predictions concerning non-linear 
density of states are made. Our results apply to other unconventional condensates such as d-wave 
superconductors and d-density waves as well. 
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I. INTRODUCTION 

Our understanding of the problem of dilute concentration of magnetic impurities in normal metals benefited a lot 
from exact solutions, and from reliable approximate methods 1 . Among the latter, the large- N expansion seems to 
describe successfully the low temperature properties of matter, N denotes the spin degeneracy, and 1/N can be used 
as an expansion parameter^^. On the other hand, recently considerable attention has been focused on the behaviour 
of magnetic impurities in pseudogap phases, where the conduction electron density of states varies as a power-law of 
energy around the Fermi energj^SiLSiSiifliii. The interest on this subject mainly arises due to the behaviour of the 
pseudogap and superconducting phase of high T c superconductors. The interplay between unconventional condensates 
and quantum magnetic impurities can reveal the nature of the underlying phase. Using large-iV technique, several 
interesting results were found concerning the transition from magnetic to non-magnetic phases, the quasiparticlc 
density of states etc. It was stated that particle-hole symmetry breaking is necessary to enter into the Kondo 
regime^, which condition is naturally satisfied with an infinite-U Anderson impurity. 

Experimentally, from scanning tunneling microscopy studies in Bi2Sr2CaCu2Q8+ j 12113 , a strong variation of the 
electron density of states was found around the impurity site: beyond the superconducting coherence peaks, new 
structures were identified indicating the presence of impurity induces bound states. Some of the experiments can 
be explained by pure potential scattering, but other features call for theoretical works to understand the effect of 
unconventional condensates on magnetic impurities^. 

This is why we have chosen to study the interaction between unconventional density waves (UDW) and magnetic 
impurities. In UDW, the gap in the quasiparticle spectrum vanishes on certain subsets of the Fermi surface, and 
its average over the Fermi surface is zeroi^i^. This causes the lack of periodic modulation of the spin and charge 
density. Such states have been proposed to describe the pseudogap phase of high T c superconductorsi§iiLi& J the low 
temperature phase of q-(BEDT-TTF)2KHg(SCN)4 19 i 20 , the antiferromagnetic phase in UR^SL, 21 : 22 and other heavy 
fermion materials^, the charge density wave in 2H-TaSe2> 2 ^ and the pseudogap phase in transition metal oxides^. 
Due to the wavevector dependence of the gap, the transition to this phase is metal to metal instead of metal to 
insulator, as in conventional density waves (with constant energy gap). 

The paper is organized as follows: in Sec. II, we study an Anderson impurity embedded in an unconventional 
density wave in the large- N limit at T = 0. As opposed to previous treatments^, we allow for a macroscopic 
occupation of the / state. From the saddle point equations, the phase diagram is determined, and the effect of 
magnetic field is discussed. In Sec. Ill, we turn to the investigation of the properties of the impurity. Its density of 
states displays the same power-law energy dependence at low energies as that of band electrons, enhances around ± A 
and a Kondo peak shows up at positive energies. The presence of these three peaks is in accord with experimental 
findings in I^Si^CaC^Os-i-,^^. The peculiarities of the conduction electron transport lifetime are explored. The 
static spin susceptibility signals the transition from the Kondo to the decoupled free moment regime. The dynamic 
spin susceptibility exhibits the usual Kondo peak, plus an additional peak coming from the divergent peak in the 
density of states of UDW. Some generalizations to non-linear density of states are made. In spite of the different 
topology of the Fermi surfaces, our results apply to other phases with power-law density of states like in d-wave 
superconductors^ or in d-density waves^t 
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II. PHASE DIAGRAM 



The Hamiltonian describing an infinite-U Anderson impurity interacting with an unconventional density wave is 
given by: 

H = J2 [^ k )( a t m a Km ~ ak_Q, ro ak-Q,m) + A(k)(a+ m a k _ Q>m + a+_ Q m a Km ) + 

k,m 

+ Ej2tifm, (1) 

m 

where aj m and a ki?Ti are, respectively, the creation and annihilation operators of an electron of momentum k and 
spin m, —s < m < s, N = 2s + 1. Similarly, /+ and / m creates and annihilates an electron on the localized E 
level, b + and 6 are the slave boson operators, responsible for the hole states^. In a sum with prime, k x runs from 
to 2kp (kF is the Fermi wavenumber), Q = (2kp , ir/b, tt/c) is the best nesting vector. A(k) = Asin(ofcj,) is the 
unconventional density wave order parameter, and is taken to be real for simplicity. A gap with cos(bk y ), or k y — > k z 
replacement would yield to the same results. Our system is based on an orthogonal lattice, with lattice constants a, b, c 
toward direction x, y, z. The system is anisotropic, the quasi-one-dimensional direction is the x axis. The a electron 
system possesses a density of states, which varies linearly with energy around the Fermi energy: N(u>)/No — \oj\/A 
for \uj\ <C A, No is the normal state density of states per spin at the Fermi energjii^. The linearized kinetic-energy 
spectrum is £(k) = vp(k x — kp), and the dispersion in the other directions can be neglected for practical purposes. 
By introducing the spinor 

the single-particle thermal Green's function of the a electrons without hybridization is obtained from Eq. as 

G(k,iuj n ) = - / dr(T T *(k,m,r)*+(k,77i,0)) ff e^ r = ^-^pg-piACk)]- 1 , (3) 
Jo 

where u n is the fermionic Matsubara frequency, piS (i = 1,2,3) are the usual Pauli matrices acting on momentum 
space. 

The Hamiltonian should be restricted to the subspace 

^tif m + b+b = Q (4) 

m 

where in the true large- -/V limit Q grows extensively with N, i.e. Q = Nq, < q < 1. In the original formulation 
of the model, Q = 1, but in order to preserve a macroscopic occupation of the / state, the above generalization was 
found to be useful^. It has been argued that the properties of the spin-i Anderson model are best represented by 
q = 1/2 in the large- N limit^. Within the mean-field approximation, the slave-boson operators are replaced by their 
expectation value, bo = (b) / y/~N, and the constraint is satisfied by introducing a Lagrange multiplier A: 



-7= ( ( a iU 



a 



k-Q,m 



)f m b + + (<z k , m + a k _Q jm ) bj 



H = [C(k)(o k " jm a k ,m ~ a k -q, m ak-Q,m) + A(k)(a£ m a k _ Q:m + a+_ Q m a k , m ) + 

k,m 

Vb o (« m + <-q, J/m + (a k ,™ + a k _ Q , m ))] +(E + X)J2 tifm + NX(b 2 - q), 

m 



(•5) 



It has been shown^ that the slave boson mean-field approximation produces the correct low energy physics of the 
conventional Anderson impurity model in the entire parameter range at T = 0. The value of A and bo is determined 
by minimizing the free energy of the system with respect to therr&S. As a result, the saddle point equations at T = 
are given by: 

,2 1 1 f° E + X 

On = a - - H — / = -ax, (6) 



7T Jo x 2 {l + ta{x)f + {E + \y 



OnA = — 



x 2 (l + Ta(x))a(x)T 
7T Jo x 2 (l + Ta(x)) 2 + (E + A) 



-dx, 



(7) 
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where a(x) = 2K(1/ yl + (x / A) 2 ) / ^/ A 2 + x 2 tt for x < W, otherwise, K{z) is the complete elliptic integral of the 
first kind. The case of a magnetic impurity embedded in a normal metallic host can be studied with a{x) = 1/\x\. 
W = vpfcF is half of the bandwidth, T = V 2 ttNq, T = Tb^. The critical value of E is obtained as 



E c 



7T \ A 



(8) 



in the weak coupling limit (W 3> A). Below this value, only the trivial solution of the saddle point equations exist, 
namely bo — and A = —E, and no Kondo effect occurs. To study the more general case, when the density of states 
varies as a power-law of energy, one has to assume A(k) = A| sin(6fc !y )|''sign(sin(6fc. y )), < r < oo. The sign function 
assures the vanishing average of the gap over the Fermi surface, and the r exponent results in a density of states as 
N(lu)/N = \Lu/A\ 1 / r T(l/2r)/^rT(l + r/2r) for \uj/A\ 1 / r < 1, and T(x) is the complete Gamma function. With 
this, the above condition is modified as 



E r = 



r [2 r+1 w 

— In 

7T V A 



(9) 



The higher the value of r, the lower the critical / level energy, and by letting r — > oo, the case of a normal metal 
can be reached with a constant density of states, when E c — > — oo, which means, that Kondo effect is always present 
at T = 0. The saddle point equations have been solved numerically, and the results are shown in Figs. ^ and [3 for 
q = 1/2 and q = 1/6. For comparison, we also show the results assuming a normal metallic host. As T decreases, bo 
increases much sharper, and reaches its maximum value ^fq rapidly. In spite of the different constraints (q = 1/2 and 
1/6), the similarity between the figures is striking, indicating that the results hardly depend on the filling factors. As 
to A, it equals to \E\ below its critical value, and approaches as the / level energy is further increased. 




E/A 

FIG. 1: The E dependence of the order parameter (the expectation value of the slave-boson operator) is plotted for q = 1/2, 
W = 20A, r/A = 0.1, 0.5, 1 and 2 from right to left (solid line), while the dashed line represents the behaviour of an Anderson 
impurity embedded in a normal metallic host. 

The effect of magnetic field can readily be incorporated into the theory by adding the Zeeman term to the Hamil- 
tonian: 



(10) 



where the appropriate choice of t m can represent the two extreme limits 2 ^: (1) the straightforward generalization of 
the SU(2) model, where t m = m/N and the magnetic field completely lifts the impurity degeneracy, or (2), when 
t m = sign(m), where two iV/2-fold degenerate levels are introduced by the field. For large fields, the critical / level 
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E/A 

FIG. 2: The E dependence of the order parameter is shown for q = 1/6, W = 20A, T/A = 0.1, 0.5, 1 and 2 from right to left 
(solid line), while the dashed line represents the behaviour of an Anderson impurity embedded in a normal metallic host. 



energy is obtained as 

q J for case (1) 

sign I - — q J for case (2) 

where _E c o is the critical / level energy without magnetic field. For q < 1/2, which is thought to be the physically 
motivated case, an applied magnetic field enhances the critical E, hence destroys the Kondo region. This follows 
naturally from the fact, that we have N/2 states below E generated by the magnetic field, so the highest occupied 
state at T = is below the actual E. On the other hand, for q > 1/2, by the use a magnetic field, the system can be 
driven back into the strong coupling regime, but the physical realization of such a situation (i.e. q > 1/2) seems to be 
doubtful. Here, the highest occupied state generated by the magnetic field is above E, which makes the hybridization 
possible again. 



III. DENSITY OF STATES, SPIN SUSCEPTIBILITY 

The Green's function of the / electrons reads as 

G f {iuj n ) 



1 



(12) 



iuj n — E — A + iuj n Ta{ui n ) 

where a(x) was defined below Eq. The self energy of the / electrons along the real frequency axis is obtained as 

2f 



+ id) = — - 



e(A-H) I^k | 



- 



+ i 



K 



6(M - A)iK 



(13) 



which exhibits the marginal Fermi liquid behaviour discovered by Zhang et al. 8 : for small frequencies, the self 
energy varies as S t (u < A) = -f(2wln(4A/M)/7r + i\u)\)/A. For r ^ 1 exponent, ReE/ ~ sign(w)M min ( 1 < 1 / r \ 
ImEj ~ M 1 / 1 ". In general, the imaginary part of the self energy is directly proportional to the a electron density of 
states. From the self energy, the real / electron density of states, which is the imaginary part of the Fourier transform 
of -(T T 6+(T)/ m (T)/+(0)6(0)>4, reads as 



Im£/(w) 



7T (w-S-A-ReS / (w)) 2 + (ImE / (w)) 2 ' 



(14) 
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which has the same low energy behaviour as the a electron density of states, namely around the Fermi energy 
Pf{u>) ~ Iw/AI 1 / 1 ". Also around ±A, it is expected to increase sharply due to the divergence of the UDW density of 
states, but the divergent peak is suppressed since the denominator in Eq. 1141 also diverges, as can be seen in Fig. [3] 
In fact, p/(±A) = 0, which is invisible in Fig. [3] due to its scale. The same phenomenon was discussed in charge 
density waves in the presence of a non-magnetic impurity 29 . Above i? c , the large, asymmetric Kondo peak appears, 
and moves to higher frequencies, even above A, as E is further increased. 
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FIG. 3: The / electron density of states is shown for q = 1/4 for E slightly above E c (solid line), E = (dashed line) and 
E — \E C \ (dashed-dotted line). 



The T-matrix of the conduction (a) electrons is obtained as 

V 2 h 2 

T{u + l 5) = — (15) 

U) — E — A — ljf(uj) 

whose poles determine the impurity bound states or resonances^. By closely examining the denominator of Eq. 1151 
we find resonances at the — *7, can be identified with the Kondo temperatur e) 2 ! 14 , as was done in similar 
treatments of the large- N limit of the Kondo mode&SiiS, and 7 is the lifetime broadening of the / level*, and is 
shown in Fig. 0] The Kondo temperature increases almost linearly with E, 7 reaches its maximum around E = A. 
For large E, the resonance appears at E — iqT, identically to the case of a normal conduction band*. Alternatively, 
one can define the Kondo energy scale 2 * 4 as \/Q 2 K + 7 2 , but it hardly differs from Qk ■ 

The imaginary part of the T-matrix, which gives the inverse transport lifetime of the conduction electrons, is 
proportional to the real / electron density of states, hence for small energies r ~ M" 1 . Its detection is a very 
subtle issue, since the number of states having this behaviour is N(u>) ~ |u>|, hence N(uj)t ~ const. The same 
cancellation is found for the r^l case. Also r is expected to be enhanced around ±A as ln(2\/2/ a/|1 — |o;/A||). The 
number of possible states diverges in the same manner as r, resulting in N(uj)t ~ In 2 - \uj/A\\. For non-linear 
density of states, similarly to the case of r = 1, both r and N(u>) diverges around the gap maximum, which gives 
N(oj)t ~ In 2 |1 — Iw/AH 1 / 2 ^. This behaviour could be checked by experiments, although its detection is not an easy 
task, because such excitation are suppressed at low temperatures. With the help of the T-matrix, the local UDW 
density of states can be determined at the impurity site: 

iVoP/H^ - E - A) 2 tt 

PaH = ^ 1 (16) 

and is shown in Fig. [S] Slightly above E c , a sharp resonance at flu modifies the density of states. By increasing E, 
the Kondo peak is broadened, and the new zero at u> = E + A can completely suppress the peak at uj = A. 
Using Eq. the low temperature specific heat is obtained as 



(17) 
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E/A 

FIG. 4: The Kondo scale, Q x is plotted for q = 1/2, V/A = 0.1, 0.5, 1 and 2 from right to left (solid line). The dashed line 
represents the behaviour of a normal metallic host. Note that Qk never vanishes in this case. The inset shows the lifetime 
broadening, 7, together with normal host results. 
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FIG. 5: The local density of states (at the impurity site) of UDW is plotted for q = 1/4 for E slightly above E c (solid line), 
E — (dashed line) and E — \E C \ (dashed-dotted line). 



which has the same temperature dependence as that of the density wave. For a non- linear density of states, the specific 
heat varies as C(T) ~ T , ( 1 + r )/ r . Also from the behaviour of p/(u>), the low temperature behaviour of the resistivity 
can be calculated^, and is obtained as R(T) — R{0) ~ — T 2 ln 2 (T). For general r / 1 exponents, it decreases as 

_jimin(2,2/r) 

The impurity spin susceptibility is calculated from 

x(ium) = —plffT/lGfjiun +iu> m )Gf(iu> n ), (18) 

n 

whose static limit at T = reads as 

,,„) . All p«-«+* _ i r .(*+*>' a , (19) 

W L// V E + X irj (x 2 (l + Ta{x)Y + (^ + A) 2 ) 2 / V ' 



which is shown in Fig. Altogether with the corresponding results in a normal metallic host. The susceptibility diverges 
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at E Cl indicating the transition from the mixed valence to the decoupled free spin region, since in the latter region, it 
behaves as X (T) = t*l ff q(l - q)/T. 
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FIG. 6: The E dependence of static spin susceptibility is shown on a semilogarithmic scale for q = 1/2, W = 20A, T/A = 0.1, 
0.5, 1 and 2 from right to left (solid line), while the dashed line accounts for the susceptibility of an Anderson impurity embedded 
in a normal metallic host. 

The imaginary part of the impurity spin susceptibility should be readily accessible by neutron scattering experi- 
ments, and is evaluated at T = as 

Imx(u) = -^4— p f (x)p f (x-u)dx. (20) 
°o Jo 

For small frequencies, it behaves as Imx(w) = T 2 lu 3 /6ir(E + A) 4 A 2 , while for non-linear density of states, it increases 
as ui 1+2 / r . Beyond the usual Kondo peak occurring around £Ik, an additional sharp peak shows up at A + f^K 
corresponding to excitations around the UDW gap. As E increases, the latter becomes dominant, but its sharpness 
is smeared, as can be seen in Fig. [7] In a pure UDW, only one single peak is expected in Imx(u;) located at 2A, so 
the presence or absence of the extra two peaks can identify the presence or absence of magnetic impurities. 

IV. CONCLUSION 

In summary we have studied the screening of a magnetic impurity in unconventional density waves. We considered 
an infinite-U Anderson impurity in the large- N limit, when the occupation of the / level increases extensively with N. 
Kondo effect only occurs, if the / level energy or the hybridization matrix element exceeds a certain value, and the 
result are almost independent of the chosen value of the / level occupation. The impurity density of states exhibits the 
UDW coherence peaks at ±A, and the usual Kondo peak. The conduction electron density of states at the impurity 
site displays similar behaviour. The transport lifetime of the conduction electrons is determined from the T-matrix, 
and is found to diverge at u) = 0, and ±A. The former is suppressed due to the small number of states possessing this 
kind of behaviour. The static impurity spin susceptibility diverges at the critical / level energy, where the transition 
from Kondo to weak-coupling takes place. The dynamic spin susceptibility is enhanced at the Kondo energy, but an 
additional sharp resonance is noticed belonging to excitations from the gap maximum of UDW states to the Kondo 
peak. These features can help to identify the magnetic or non- magnetic nature of impurities in pseudogap phases. 
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FIG. 7: The evolution of the dynamic spin susceptibility for various values of E (from E — — 1.25A by 0.1A steps from left 
to right) is shown for q = 1/6, W = 20A, T/A = 1. Two distinct peaks show up above E c , and merge as one passes into the 
Kondo regime with increasing E. The circles denote E £s 0. 
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